Rate dependent shear bands in a shear transformation zone model of amorphous sohds 
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We use Shear Transformation Zone (STZ) theory to develop a deformation map for amorphous 
solids as a function of the imposed shear rate and initial material preparation. The STZ formulation 
incorporates recent simulation results [Haxton and Liu, PRL 99 195701 (2007)] showing that the 
steady state effective temperature is rate dependent. The resulting model predicts a wide range of 
deformation behavior as a function of the initial conditions, including homogeneous deformation, 
broad shear bands, extremely thin shear bands, and the onset of material failure. In particular, the 
STZ model predicts homogeneous deformation for shorter quench times and lower strain rates, and 
inhomogeneous deformation for longer quench times and higher strain rates. The location of the 
transition between homogeneous and inhomogeneous flow on the deformation map is determined in 
part by the steady state effective temperature, which is likely material dependent. This model also 
suggests that material failure occurs due to a runaway feedback between shear heating and the local 
disorder, and provides an explanation for the thickness of shear bands near the onset of material 
failure. We flnd that this model, which resolves dynamics within a sheared material interface, 
predicts that the stress weakens with strain much more rapidly than a similar model which uses a 
single state variable to specify internal dynamics on the interface. 

PACS numbers: 83.60.Rs,83.50.-v,46.35.+z,81.40.Lm, 62.20.F 
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I. INTRODUCTION 



Amorphous solids such as foams, dense colloids, bulk 
metallic glasses, and granular fault gouge are ubiquitous 
in engineering applications and natural systems. Al- 
though these materials exhibit a yield stress on exper- 
imental time scales, they flow, deform, and fail in a man- 
ner which is different from crystalline solids or Newtonian 
fluids. Many of these materials undergo strain localiza- 
tion, where a small region deforms much more rapidly 
than adjacent regions. For example, bulk metallic glasses 
develop very thin shear bands [l|, Q , fault gouge in earth- 
quake faults develops a prominent fracture surface that 
accommodates most of the slip jSj] , and colloidal systems 
develop broad shear bands [1]. 

Surprisingly, the mechanisms that lead to this strain 
localization have remained elusive. An early theory of 
shear banding [Bj suggests that a small increase in the 
thermal temperature lowers the viscosity, resulting in 
more rapid deformation and a local increase in temper- 
ature. However, Lewandowski and Greer showed that 
shear bands in bulk metallic glasses can not be explained 
by adiabatic thermal effects [6|]. Although thermal heat- 
ing must play a role at high strain rates, it does not 
appear to govern the formation of shear bands in many 
materials. 

In a recent paper Q , we found that at low strain rates 
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the Shear Transformation Zone (STZ) theory for amor- 
phous solids predicts shear band formation. These bands 
are generated by feedback between the local strain rate 
and the configurational disorder of an amorphous pack- 
ing. In a separate paper Q, we used STZ theory to fit 
data over a wide range of strain rates from a simulation of 
glassy disks by Haxton and Liu @ . These data indicate 
that dynamics of the glassy material change dramatically 
at large strain rates. 

Experimental studies of bulk metallic glasses driven at 
a wide range of strain rates also show that deformation is 
strain rate dependent. Homogeneous deformation is seen 
at low strain rates and inhomogeneous flows dominate at 
large strain rates [l3, [IH ■ 

In this paper, we study deformation in the STZ model 
as a function of the initial material preparation and ex- 
ternally imposed strain rate, and include the rate depen- 
dence observed by Haxton and Liu in our model. This 
STZ formulation is valid for a large range of strain rates 
and predicts four different types of deformation behavior: 
homogeneous deformation, thick disorder limited shear 
bands, thin diffusion limited shear bands, and material 
failure. We will discuss these various types of deforma- 
tion in detail in the following sections. We describe the 
mechanism that generates shear bands and the processes 
that determine the thickness and longevity of these in- 
homogeneous flows. We numerically integrate the STZ 
equations to produce a deformation map for the model 
glass simulated by Haxton and Liu that shows the type 
of deformation or failure predicted as a function of the 
initial conditions. 

Shear bands in this STZ formulation occur due to a 
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feedback between the effective temperature, which de- 
scribes the configurational disorder in a glassy or jammed 
material, and the local strain rate. In a sheared, steady 
state, non equilibrium amorphous material, the effective 
temperature can be calculated by measuring the fluctu- 
ations and linear response of an observable such as the 
pressure and applying the fluctuation-dissipation theo- 
rem (FDT) ,9, 12. 13]. Ono, et al. have shown that in 
several simulated foams, measurements of different ob- 
servables yield a single, rate dependent steady state effec- 
tive temperature which is distinct from the thermal tem- 
perature [l^. In addition, these authors show that the 
FDT effective temperature is consistent with an entropic 
definition: the effective temperature is the derivative of 
the configurational entropy with respect to the (poten- 
tial) energy. This definition suggests that slow steady 
shearing causes the material to ergodically explore all 
possible configurational packings, and therefore the sys- 
tem maximizes a configurational entropy. 

The STZ model describes plastic deformation in amor- 
phous material as a series of localized plastic events that 
occur in susceptible regions, or zones [ll,[l^,[l^. Follow- 
ing Falk and Langer, we model STZs as bistable regions 
that are more likely than the surrounding material to 
deform under stress, and are created and annihilated as 
energy is dissipated in the system [T^. This model has 
successfully been used to describe bulk metallic glasses, 
thin filrns, and hard spheres in several different geome- 
tries 0, Ha, m, Hill. 

An important feature of all STZ formulations is that 
the zones are activated by an effective temperature or 
free volume, and there is a feedback between packing 
structure and deformation. In particular, we postulate 
that an STZ is an unlikely, high energy configuration of 
an amorphous packing. Because the effective temper- 
ature governs the statistics of configurational packings, 
the steady state density of STZs should correspond to a 
Boltzmann factor. 
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where A is proportional to the steady state STZ density, 
Ez is an activation energy, fcs is the Boltzmann constant, 
Te// is the effective temperature, and x = Ez/{kBTeff) 
is a dimensionless effective temperature [l7| . 

Because plastic deformation occurs only at these STZs, 
the plastic strain rate in simple shear is proportional to 
this density, and in many situations can be written as 
follows: 
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where 7 is strain, /(s) is a function of the deviatoric 
stress s and tq is an internal timescale such as the phonon 
frequency. For the remainder of this paper we will refer 
to the dimensionless plastic strain rate, q = 'jpiTo and 
q = (Vq/L)tq is the imposed average strain rate times 
the STZ time scale. 



For completeness. Appendix [XI reviews the STZ model 
and defines the exact equations and parameters used in 
this work. The model is vastly simplified by focusing on 
materials in a simple shear geometry at thermal temper- 
atures far below the glass transition temperature, and 
can be summarized by two equations. 

The first specifies that elastic deformation increases 
the shear stress, while plastic deformation decreases it: 
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= M* 
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where 11^ is the ratio of the elastic modulus to the yield 
stress, and A is the spatial average of the STZ density 
A = exp(— 1/x). In writing this equation we have as- 
sumed that the deviatoric stress is constant across the 
width of the sample, as discussed in the Appendix. 

A second equation specifies that the effective temper- 
ature approaches its steady state value, xil) as plastic 
work is dissipated, and it also diffuses: 

«7 CQSoq XW oy^ 

where cq is a dimensionless specific heat, a is a diffusion 
length scale on the order of an STZ radius, and sq is the 
stress threshold for the onset of plastic deformation. For 
the remainder of this paper, a symbol with an overline 
denotes a quantity that is a spatial average or constant 
as a function of position. The function /(s) is described 
in the Appendix. 

The remaining unspecified component of Eq. [4] is the 
steady state effective temperature x- This material de- 
pendent parameter captures the physically intuitive idea 
that there is a maximum possible disorder attained in 
a sheared amorphous packing; above x no heat can be 
dissipated in the configurational degrees of freedom. Re- 
cent simulations by Haxton and Liu Q suggest that the 
steady state effective temperature is dependent on the 
plastic strain rate, denoted q. In steady state, the in- 
verse strain rate 1/q can be viewed as a function of the 
steady state effective temperature x, much like the vis- 
cosity is a function of the thermal temperature T. We 
have shown Q that the the glassy steady state effective 
temperature is well fit by the following functional form: 
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A discussion of this effective temperature glass transition 
is given in Q, we adopt the super- Arrhenius function 
aeff identified in that paper: 
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where xo is the thermal glass transition temperature, XA 
is the temperature at which the system begins to exhibit 
Arrhenius behavior, and b, xi £^re fit parameters. The 
super- Arrhenius component (Eq. ^) ensures that the 
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effective temperature approaches a constant, %0j as the 
strain rate approaches zero, which is seen in simulations. 

As noted above, this paper focuses on deformation in 
geometries where the equiUbrium shear stress is constant 
across the sample. Although there are many interest- 
ing geometries with shear stress gradients (and STZ the- 
ory can explain these shear bands ^2d\), often in these 
cases stress effects compete with internal structural ef- 
fects and complicate the analysis. In contrast, exper- 
iments in simple shear geometries exhibit strain local- 
ization even though the equilibrium stress is constant in 
space, indicating that some property of the internal state 
governs shear band formation. 

Shear banding is caused by the coupling in Eq. ^ be- 
tween a configurational or structural parameter and the 
plastic rate of deformation. Although the stress equili- 
brates quickly, the effective temperature dynamics con- 
tinue to evolve over much longer timescales. Materials 
that develop long-lived shear bands exhibit a very differ- 
ent macroscopic rheology compared to those that deform 
homogeneously [2lj. The type of deformation has im- 
portant implications for the macroscopic material and 
frictional properties — systems that localize also weaken 
rapidly. 

The remainder of this paper is organized as follows. 
In Section |TT] we study the stability of the model given 
by Eqs. ^ and ^ with respect to perturbations. We 
numerically integrate the STZ equations to validate our 
analytic stability results in Section Hill and study the dif- 
ferent types of deformation that persist for long times in 
Section [iVl Section IVl concludes with a discussion of our 
results and open questions. 



II. STABILITY ANALYSIS 

We now study shear band formation for systems with a 
rate dependent steady state effective temperature xil)- 
The fact that shear bands persist in the STZ model is 
somewhat surprising, given that the only stationary state 
of Eqs. ^ and (|3|) is homogeneous deformation. In the 
following sections, we provide an explanation for shear 
band formation and evolution across a wide range of 
strain rates. 

We emphasize that all of the shear bands discussed be- 
low are transient phenomena ~ at very large strains the 
shear bands diffuse across the entire width of the mate- 
rial and deformation once again becomes homogeneous. 
However, because the internal dynamics of the configu- 
rational degrees of freedom can be very slow compared 
to the stress evolution timescale, shear bands persist for 
very long times. For example, the shear stress in numer- 
ical solutions to the STZ equations generally appears to 
reach a steady state within a few percent strain but the 
effective temperature profile often remains inhomogenous 
for more than 20% strain. 

First, the evolution equations for x E^nd s can be writ- 



ten as follows: 
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where the (•) operator indicates a derivative with re- 
spect to time in units of strain. The function xil) — 
x(2j{s)e~^l^) depends implicitly on the stress and the 
effective temperature. Then the Jacobian J is given by: 
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The term W{6x) is a spatial integral over one period of 
the perturbation function; it selects only the zero wave 
number component of the perturbing function because 
the other components must satisfy periodic boundary 
conditions: 
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1 for fc = 0, 
for k^Q. 
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Details can be found in [7| . We will use the term "zero- 
mean perturbation" to refer to a perturbing function with 
a vanishing fc = component that does not change the 
average value of the underlying function across the width 
of the sample. 

For a specified externally imposed strain rate q — 
toVq/L, the steady state solution to Eqs. [7] and [8] is 
X = x{q): and Sss is given by the solution to the alge- 
braic equation 1 = 2/(sss) exp[— l/x(g)]/g. The Jaco- 
bian evaluated at this steady state is greatly simplified: 
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The upper right entry ( J12) vanishes for zero-mean per- 
turbations, so Jacobian is lower triangular, and the di- 
agonal values alone specify stability. Jn is strictly nega- 
tive because / is a monotonically increasing function of 
s, and therefore the system is linearly stable with respect 
to zero-mean perturbations from a homogeneous steady 
state if J22 is negative, or equivalently, dx/dx < 1- 

The partial derivative dx/dx can be written as follows: 



dx 
dx 



x'{q)q 



(19) 



where x' il) is determined by the derivative of Eq. [5l 

This clarifies the physical meaning of the criterion 
dx/dx < 1. At low strain rates the super- Arrhenius 
function aeff ensures that x changes slowly with q and 
the partial derivative in Eq. [19] is significantly smaller 
than unity. In [3] we used the approximation that xio) 
is constant at low strain rates; in this case dx/dx is triv- 
ially zero and the STZ equations are always stable with 
respect to perturbations in steady state. 

At high strain rates, where aeff is zero and x{q) = 
A/\n{qo/q), it can be shown that dx/dx — 1/A. In this 
regime, the equations are stable with respect to pertur- 
bations in steady state if and only if the normalized ac- 
tivation energy A is greater than unity. 

The material dependent parameter A is greater than 
unity in all homogeneously deforming materials where 
the steady state stress increases with increasing imposed 
strain rate. To see this, note that Eq.[5]can be rewritten 
as follows: 



ln(g/(7*)' 



(20) 



where A^ is an activation energy and is a dimensionless 
constant. In steady state, x — X everywhere, so inserting 
Eq. [50] into Eq. [2] and taking the derivative with respect 
to q results in the following equation for s'(g), which is 
the derivative of the steady state stress with respect to 
the imposed driving rate: 



s'{q) = [-l + A4q)+q\n{q/q,)A:{q))]x 



-1/A,{q) 



(21) 



The function / increases monotonically with s. There- 
fore s' (q) is positive and the steady state stress increases 
with increasing imposed strain rate whenever the first 
factor in Eq. [21] is positive. In high strain regimes, the 
activation energy is a constant, A^{q) = A, and the ma- 
terial is steady state strengthening whenever A > 1. 

In most simulations of disordered disks with hard 
core repulsion (including simulation data by Haxton and 
Liu 0] and Shi et al. [2l|), the steady state stress in- 
creases with increasing strain rate, so that dx/dx is less 
than unity and the system is stable with respect to zero 
mean perturbations. Materials with a steady state stress 
that increases as a function of driving rate are called 



rate strengthening. In contrast, several experiments that 
study friction of granular fault gouge find rate weakening 
behavior; the steady state stress decreases as the imposed 
strain rate increases ^22!]. These materials might be sus- 
ceptible to shear banding even in steady state. In this 
paper we focus on shear banding in materials that are 
rate strengthening, but rate weakening materials are an 
interesting avenue of further research. 

STZ theory predicts that rate strengthening materi- 
als should continue to deform homogeneously once they 
attain steady state. However, for certain initial con- 
ditions simulated rate strengthening materials develop 
shear bands before they reach this homogeneous steady 
state. This leads us to study the stability of perturba- 
tions at each point in time along a time-varying trajec- 
tory, assuming that the system deforms homogeneously 
up to that point. This is a "frozen time" stability anal- 
ysis. 

Due to the integral in Eq. ([7]) , spatial perturbations to 
X with zero mean do not change the equation for s to first 
order in (5x : s{s,x + ^x) = K^iX) + ^{^X^)- In this case 
the linear stability is determined solely by the equation 
for X and the sign of J22- Above the yield stress, x and 
/(s) are strictly greater than zero and therefore J22 is 
negative (and the trajectory is linearly stable) whenever 
the following criterion is met: 



x>Xcrit = -^{~x + x^ + qx'{q)+ 

^ix' + i-X + r+qx'iqW) ■ (22) 



By inserting values for x{q) into Eq. (|22[) . we derive a 
linear stability prediction for the boundary between these 
two regimes as a function of the average initial effective 
temperature x and the applied strain rate o. We choose 
xiq) to fit the data from Haxton and Liu [1, 

Unlike linear stability analysis for steady states, frozen 
time stability analysis for time varying trajectories does 
not predict the final state of the system. It provides an 
indication that a transient instability is possible, but it 
does not specify global stability. The frozen time analy- 
sis is accurate only when the perturbations grow rapidly 
compared to the growth of the underlying trajectory. 
We therefore use the more general localization ratio TZ 
to characterize the transient instability. First discussed 
in 0, this ratio compares the growth rate of pertur- 
bations (determined by frozen-time stability analysis) to 
the growth of the underlying trajectory: 



n = Sx 



exp[J22(s )Ai/2]J22(Sm,X) 
X{Sm,x) 



(23) 



where 5x is the magnitude of the perturbation, Xini is the 
initial effective temperature, Sm is the approximate max- 
imum shear stress given by the solution of the equation 
q = 2f{sm) exp[— l/xi„i], and At is the approximate 
time in units of strain it takes to achieve the stress maxi- 
mum. Localization occurs when the rate at which heat is 
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FIG. 1: (color online) Predicted deformation map based on 
initial conditions only. The solid line marked with circles 
(black) represents the frozen time linear stability criterion 
predicted by Eq. ((22}, which does not take into account fi- 
nite amplitude perturbations. 7?. is a more accurate general- 
ized stability criterion that takes into account finite amplitude 
perturbations. A line with constant localization ratio TZ = 1 
is marked with crosses(blue). Above the line, TZ < 1 and 
homogeneous deformation is predicted, while below the line 
TZ > 1 and strain localization is predicted. For reference, the 
upper dashed line shows the xil) fit to the data from Haxton 
and Liu 9]. 



dissipated inside the band is larger than the rate outside 
the band; in this case TZ is greater than unity. 

The localization ratio given by Eq. depends on the 
magnitude of the perturbation Sx- In the low strain rate 
limit, we have systematically studied the localization ra- 
tio 7?. as a function of perturbation amplitudes and found 
that it accurately predicts that larger perturbations lead 
to enhanced localization ^TJ . For simplicity, we have cho- 
sen Sx to be 5 % of the average value of the effective 
temperature, which is consistent with perturbations to 
the potential energy per atom for a Lennard Jones glass 
calculated by Shi, et al. [2l[. Systematically studying the 
effects of perturbation magnitude as a function of strain 
rate is beyond the scope of this paper. 

Figure [T] is a deformation map that predicts the type 
of flow as a function of the initial conditions for the sim- 
ulated glassy material studied by Haxton and Liu 0. 
The bold line is the linear stability criterion defined by 
Eq. ^F^. Because a frozen time analysis does not take 
into account finite amplitude perturbations or the growth 
rate of the underlying trajectory, we use the localiza- 
tion ratio, TZ, to predict localization. Using Eq. E51 
we calculate TZ for each set of initial conditions with 
5x = 0.05 X x^m, and M = 0.06= 6 % strain. The 
line marked with crosses in Fig. [1] corresponds to a line 
with constant TZ — 1. Localization is expected below this 
line, where TZ > 1, and homogeneous flow above it. 



III. NUMERICAL SOLUTIONS TO STZ 
EQUATIONS 

To check these analytic predictions, we numerically in- 
tegrate the STZ partial differential equations. The nu- 
merical solutions exhibit three broad categories of de- 
formation behavior: homogeneous deformation, shear 
bands, and melting or failure. This section discusses 
qualitative features of each kind of deformation, while 
Section HVl develops a deformation map using a quantita- 
tive criterion for each category and discusses macroscopic 
implications. 

To resolve extreme localization, we use an irregular 
mesh and a combination of fixed-step and adaptive-step 
finite difference methods. For each pair of initial condi- 
tions, the average initial effective temperature Xim and 
the externally applied strain rate q = ro(Vb/L), we nu- 
merically integrate the STZ equations (Eqs. ([3|) and (j4|)) 
from to 20 % strain. The initial effective temperature 
function Xiiii (y) is a constant perturbed by a hyperbolic 
secant function of height Sx and width L/10, normalized 
so that its average is Xim, and the initial shear stress 
is 0.0001. All stresses are in units of the yield stress Sy 
unless otherwise noted. 

For comparison, we also numerically integrate a single 
degree of freedom STZ model, where the effective tem- 
perature is constrained to be constant across the width 
of the material, and no perturbations are permitted. The 
system of ordinary differential equations given by Eqs. ([3|) 
and ([4]) (with no diffusion) is integrated numerically in 
time using the same average initial conditions as the STZ 
PDF. 

The simple ODE model cannot localize and has been 
used to describe macroscopic frictional behavior for 
boundary lubrication in thin films [l^ and on earthquake 
faults [23] . The ODE model is an example of a "rate and 
state" friction law. These laws are frequently used in geo- 
physical modeling of earthquake ruptures, and describe 
the response of a sheared frictional interface as a func- 
tion of the slip rate (or strain rate) and a single state 
variable. While the STZ PDF resolves internal dynam- 
ics of the effective temperature within the interface, the 
STZ ODE is constant across the interface. Comparing 
the two models allows us to study the effect of small scale 
dynamics such as strain localization on model predictions 
for macroscopic behavior. 

In some simulations of the full PDF model, the steady 
state effective temperature x approaches infinity. Al- 
though the STZ model given by Eqs. ^ and ^ is stiU 
well-defined in this limit, the shear heating term becomes 
considerably amplified, indicating a situation where the 
amorphous packing becomes more and more disordered 
inside the band. In every instance where x ~^ oo, the 
shear band becomes so thin that the numerical integra- 
tion routine fails. 

We suggest that this numerical failure corresponds to 
material failure. The smallest length scale in the model 
is a, the diffusion length scale which is on the order of 
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the radius of an STZ. We do not expect the STZ model 
to hold at length scales smaller than a, and because our 
numerical mesh is fine enough to resolve a band ten times 
smaller than a, numerical failure corresponds to a shear 
band that rapidly becomes so thin that the model itself 
breaks down. 

Although the simple STZ model developed here does 
not specify the rheology at strain rates above this "melt- 
ing" point, it does suggest that the sohd-hke STZ theory 
must be replaced by a liquid-like theory (such as mode 
coupling or Bagnold scaling) inside these bands. There- 
fore, integration of the STZ model indicates that when 
the disorder temperature approaches infinity, the mate- 
rial can no longer support a static shear stress; it liquefies 
and fails. 

In simulations where the effective temperature remains 
finite, we numerically track the shear stress s and the 
effective temperature field xiv) functions of time, or 
equivalently, strain. In each case, the stress first responds 
elastically, and then begins to deform plastically above 
the yield stress Sy. As plastic deformation increases the 
effective temperature, the material softens and the stress 
relaxes to its flowing value, s/. The dashed blue line in 
Fig. [2] is a plot of the stress vs. the strain for a numerical 
solution to the STZ PDE model. Initial conditions are 
such that the material is highly unstable with respect to 
shear bands. For comparison, the dash-dotted (magenta) 
line in Fig. [5] shows the solution to the STZ ODE model, 
which is a rate and state law with a single internal state 
variable. The STZ PDE solution develops a shear band 
and weakens much more rapidly than a rate and state 
model with similar initial conditions. 



1.5 



S3 

N 



0.5 



o 



- K 





ODE 




- PDE 


o 


0.0239 


< 


0.0272 


□ 


0.0281 




0.0289 




0.0322 




0.0447 


+ 


0.0781 



0.02 



0.04 0.06 
strain 



0.08 



0.1 



FIG. 2: (color online) Shear stress s vs. strain calculated 
by numerically integrating the STZ equations of motion with 
initial conditions Xi-m ~ 0.0674, and imposed strain rate 
q = 1.015 X IQ-^ The dashed (blue) curve represents the 
solution to the perturbed STZ PDE model, while the dash- 
dotted (magenta) curve represents a solution to the STZ ODE 
model with the same average initial conditions. The colored 
symbols correspond to the plots shown in Figs. |3] and (5) At 
about 2% strain, the perturbed system begins to localize and 
weakens much more rapidly than the magenta curve. 



0.22 
0.2 
2 0.18 

B 

S 0.16 
u 



g 0.14 

1 0.12 
o 

2 0.1 

0.08 



0.06 



(a) 





0.0000 




0.0239 




0.0272 




0.0281 




0.0289 




0.0322 




0.0447 


+ 


0.0781 




0.1031 





-1 



-0.5 



0.5 



y - position 



FIG. 3: Diffusion-limited shear band(color online) (a) 
Time series of the effective temperature as a function of posi- 
tion for a material with initial conditions Xini = 0.0674, and 
imposed strain rate q = 1. 01510"®. Each line represents the 
effective temperature field as a function of position at a differ- 
ent time, as indicated by the legend (all times are in units of 
strain). The effective temperature field is initially a constant 
with a small perturbation centered in the middle(blue). This 
perturbation grows rapidly (green) and forms a shear band, 
which then diffuses outward slowly (red), (b) Inset shows 
the same data on a different scale. The strain associated with 
each line is also indicated in Fig. [5] localization coincides with 
rapid dynamic weakening of the shear stress. 
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FIG. 4: Average plastic strain rate, / dy ^(y)To/q as a func- 
tion of strain for the same integration data shown in Fig. [2] 
Initially the average plastic strain rate is zero during the ma- 
terial elastic response. The average plastic strain rate then 
rises rapidly during the stress overshoot, when the system re- 
leases stored elastic energy. Finally, it relaxes back to unity 
in the flowing regime, when all the strain must be accommo- 
dated plastically. 

During this initial stress response, the effective temper- 
ature field also evolves in space and time. The effective 
temperature field is initially constant with a small, cen- 
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FIG. 5: Diffusion-limited shear band(color online) (a) 
Time series of tlie normalized plastic strain rate "/ro/q as a 
function of position. The plastic strain rate is initially zero 
(during elastic deformation) but rises at the onset of plastic 
deformation and becomes very sharply peaked (green). This 
shear band is extremely narrow with a thickness of about 
0.015, which is approximately the same as the diffusion length 
scale a. As the stress relaxes the strain rate also relaxes, 
and the shear band becomes wider and less sharply peaked 
(red).(b)Inset shows magnified position and strain rate axes. 
Although the maximum strain rate in the band decays signif- 
icantly with time, it remains large (> 25 times the imposed 
strain rate). 



tered perturbation, and the field remains static during 
the elastic response. At the onset of plastic deformation, 
the effective temperature begins to rise. In systems that 
deform homogeneously the average value of the effective 
temperature rises and the perturbation dissipates, while 
in systems which develop shear bands the effective tem- 
perature rises rapidly inside the band and attains a slowly 
evolving state where the shear bands diffuse outwards. 

Figures [^a) and (b) are a series of plots of the effec- 
tive temperature as a function of position for a material 
that develops a thin shear band (later we will identify the 
initial transient as diffusion limited). Each colored line 
represents a different time in units of strain. A small ini- 
tial perturbation to the effective temperature is driven by 
a dynamic instability to a much higher value, saturating 
at X ~ 0.22, and the band then slowly diffuses outward. 

The plastic strain rate also evolves during the initial 
transient response. We first focus on the average plas- 
tic strain rate, shown in Fig. |4] as a function of time. 
At early times when the stress is below the yield stress 
Sy, the system deforms elastically and the plastic strain 
rate is zero everywhere. At the onset of plastic deforma- 
tion the average plastic strain rate increases continuously 
from zero, attains a maximum, and then relaxes back to 
the externally imposed strain rate (in the flowing regime 
all the deformation is plastic.) While the plastic strain 
rate is greater than unity, stored elastic energy is being 
dissipated. 

Although the stress is constant across the width of the 
material, regions with a higher effective temperature de- 



form more rapidly. To effectively compare strain local- 
ization at various strain rates, we plot the dimension- 
less strain rate, which is the strain rate at each location 
divided by the externally imposed strain rate. In fig- 
ures [5ja) and (b), each colored line represents a different 
time in units of strain; plots show the plastic strain rate 
as a function of position. Initially the stress is below the 
yield stress and the plastic strain rate is zero (blue line). 
Localization of the effective temperature field results in a 
very narrow peak in the strain rate field (green line). The 
strain rate in the center of the shear band is nearly 3000 
times larger than the externally imposed strain rate. As 
the stress continues to relax, the strain rate becomes less 
sharply peaked (red line). The inset plots magnify the 
position axis. 

Comparing numerical results to analytic predictions re- 
quires a method for measuring the degree of localization 
in a given numerical simulation. The degree of localiza- 
tion can be quantified using the Gini coefficient (j) [23 |. 
defined as: 



0(t) = 



1 



5]5]|P'''(y.„i)-I?^'(2/„i)|, (24) 



where {yi} is a uniform grid of n points in position space. 
The Gini coefficient is equal to zero if the material de- 
forms homogeneously and increases as the plastic strain 
rate field becomes more sharply peaked. A delta func- 
tion has a Gini coefficient of 1. During a given numerical 
simulation, the Gini coefficient 4>{t) starts out as a very 
small number and then increases rapidly as the shear 
band forms. Then, as the shear band diffuses the Gini 
coefficient decreases. Because we are focusing on the ini- 
tial transient, we first study the maximum value of the 
Gini coefficient attained during a given numerical simu- 
lation. 

Figure [nja) is an intensity plot of the maximum value 
of the Gini coefficient as a function of the average initial 
effective temperature, Xim and the natural logarithm of 
the dimensionless imposed strain rate \og{q). This defor- 
mation map indicates that material deformation grad- 
ually changes from homogeneous flow to shear banding 
as a function of the initial conditions. In figure [SJa), 
black boxes indicate that x approached infinity during a 
particular numerical integration, and the STZ solid-like 
description breaks down. 

While the Gini coefficient is a direct indicator of local- 
ization, it is perhaps a less familiar metric. For compari- 
son. Figure inUb) shows that maximum plastic strain rate 
attained in the band as a function of the initial condi- 
tions. A larger plastic strain rate is attained in a thinner, 
more localized band, and therefore Fig. is very sim- 
ilar to Figinja). Again, black boxes correspond to shear 
bands where the plastic strain rate reaches the melting 
point and the model breaks down. 
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FIG. 6; (color online) A diagram showing the degree of lo- 
calization found in by numerically integrating the STZ equa- 
tions, (a) The maximum Gini coefficient (colorbar), Eq. (I24|l . 
as a function of the average initial effective temperature and 
the externally imposed strain rate. A higher Gini coefficient 
indicates more localization, (b) The log of the maximum plas- 
tic strain rate attained in the band(colorbar) divided by the 
externally imposed strain rate. In both figures, black boxes 
correspond to numerical simulations where the magnitude of 
the strain rate was so large that x ^ cx, as discussed in the 
text. The solid white line corresponds to the predicted lo- 
calization ratio TZ = 1: Localization is expected for initial 
conditions below this line. See, e.g. results for Xini = 0.0674, 
q = 1.015 X 10"'' in Fig.4. 



A. Shear band thickness 

As mentioned earlier, the stress appears to achieve 
a steady state quickly - after less than 7 % strain all 
numerical STZ solutions have acheived a steady stress 
that changes by less than 5 % over the course of the 
remaining simulation (200 % strain). In comparison, 
the effective temperature field often remains highly local- 
ized for t > 20% strain, and broadens over much longer 
timescales than the stress. 

The goal of this section is to calculate shear band 
thickness as a function of initial conditions within the 
STZ model, and determine what sets the thickness of the 
shear bands in this model. Because localized strain states 
are transient - the bands diffuse outward over time - we 
study the model predictions for how shear band thickness 
evolves over large (20 %) strains. Importantly, many ini- 
tial conditions lead to numerical simulations and exper- 
imental materials that fail before reaching large strains. 
The following analysis identifies these events as well. 

We calculate the shear band thickness for each numer- 
ical solution at two times: the time tqmax at which the 
strain rate in the shear band attains its peak and the 
shear band thickness is minimized, and at a later time 
t = 20% strain where the stress appears to be in steady 
state. For systematic study, we specify initial conditions 
that generate only a single shear band - multiple shear 
bands are often found experimentally at higher strain 
rates and will be a topic of future study. 

We first study the shear bands at a time when the plas- 
tic strain rate is most highly localized. Let qmaxijj) be 
the normalized plastic strain rate 'j(jj,t)TQ /q evaluated 
at the time tqmax when the strain rate achieves its abso- 
lute maximum. Sqmax is the shear stress at tqmax- The 
thickness Wqmax of the shear band in a numerical solu- 
tion is defined to be the fraction of the real line between 
— 1 and 1 where the function qmaxiu) is sufficiently large: 



^qmax — / ^V- 
II 



The region / is defined as follows: 

/ = {y e [-1, 1] I qraax{y) > 1 -f 



h sup qr, 

V 



(25) 



Ay)}. (26) 



IV. DEFORMATION MAP AND 
MACROSCOPIC IMPLICATIONS 



Numerical solutions presented in the previous section 
show that transient dynamics can lead to inhomogeneous 
flows. We would like to understand how to characterize 
these flows. What type of deformation occurs as a func- 
tion of the initial conditions? If shear bands form, what 
sets their thickness? What are the implications of inho- 
mogeneous flows for macroscopic system response? 



where h is an arbitrary fraction. Although we choose 
h = l/IO, in most cases the calculated thickness is insen- 
sitive to the value of h because the strain rate function 
is sharply peaked. When the system deforms homoge- 
neously, the thickness Wqmax is not well-defined. In this 
case Eq. ((25)) becomes extremely sensitive to the fraction 
h and is no longer accurate. 

Figure [7] is a plot of the shear band thickness at the 
time of maximum strain rate Wqmax as a function of the 
initial conditions. White boxes in Figure [7] correspond 
to solutions where the maximum Gini coefficient is less 
than 0.35. These homogeneously deforming solutions do 
not have a well-defined value Wnmm-- 



9 




-13 -12 -11 -10 -9 
ln(imposed strain rate) 



& 

E 0.18 



■B 0.16 




-12 -11 -10 -9 
ln{imposcd strain rate) 



-0.25 



FIG. 7: (color online) Shading indicates the shear band thick- 
ness at the time of maximum strain rate Wqmax, Eq. (|29|l . for 
numerical STZ solutions as a function of the initial conditions 
Xini and log(g). The black boxes correspond to initial condi- 
tions for which x — + oo during an integration, while the white 
boxes correspond to initial conditions for which the flow is 
homogeneous (the maximum Gini coefficient, (j>max < 0.35) 
The color scale is set such that the maximum thickness is 0.2. 



The discussion in the previous paragraphs analyzes 
shear bands at their peak, when the plastic strain rate 
is maximized. This generally occurs at less than 7 % 
strain. Fig. [5] shows the shear band thickness at at 20 % 
strain. In each case the shear bands have become wider, 
as expected. (Note that the maximum thickness shown 
in this plot is 0.3, as compared to 0.2 in Fig. \7\) Al- 
though these systems do not acheive a stationary state, 
the slowly evolving shear band thickness is observable, 
and has been seen in molecular dynamics simulations 
where periodic boundary conditions allow the system to 
be studied at very large strains. 

Perhaps the most interesting feature of Figs. [7] and [H] 
is that the system exhibits no obviously preferred length 
scale - the shear band thickness varies continuously from 
about 0{a) ~ 0.015 to 0(1) (homogeneous flow). More- 
over, the shear band thickness increases with time. Both 
of these observations are a consequence of the fact that 
localized states are transient solutions to the equations 
of motion rather than steady state solutions. 

Although there is no preferred shear band length scale, 
the shear band thickness is reproducible; the STZ model 
generates shear bands of the same thickness given the 
same average initial conditions, even if the perturbations 
are random [T]] . In addition, simulations of Lennard Jones 
glasses generate reproducible shear band thicknesses as 
a function of time [2l| and experiments on bulk metal- 
lic glasses find a characteristic shear band thickness [ll[ . 
These results suggest that on a given observational time 
scale, the system does pick out a specific shear band 
thickness. Because it is observable and reproducible, the 
shear band thickness must evolve very slowly compared 
to the stress relaxation time scale. We exploit this fea- 



FIG. 8: (color online) Shading indicates shear band thick- 
ness at 20 % strain for numerical STZ solutions as a function 
of the initial conditions Xini a^nd log(g). The black boxes 
correspond to initial conditions for which x ~^ during an 
integration, while the white boxes correspond to initial condi- 
tions for which the flow is homogeneous (the Gini coefficient 
at 20 % strain, <p{t = 0.2) < 0.35) Note that the scale for 
this plot is larger than that in Figure 7 - the shear bands are 
significantly wider at 20 % strain than at tqmax- 

ture, showing that the STZ model singles out three dif- 
ferent deformation profiles that evolve slowly in time and 
which should therefore describe observable deformation 
modes. In addition, we discuss another state - material 
failure - where a parameter in the STZ model diverges 
and the model fails. 

We analyze Eq. [8] to determine what deformation pro- 
files generate the smallest change in the effective temper- 
ature. Because the stress is nearly stationary, flows with 
the smallest average values for x are the longest-lived 
transients and are easily observable. We show deforma- 
tion profiles for these states and develop a deformation 
map at tqmax and t — 20% strain. This provides an ex- 
planation for observed shear band thicknesses. 



B. Relaxation towards homogeneous deformation 

The first and simplest state minimizes x everywhere 
and is "homogeneous deformation." The effective tem- 
perature field is constant everywhere and equal to x{q)i 
where q is the externally imposed dimensionless strain 
rate. Since both the shear heating and diffusion terms 
are zero in Eq. ([5]), this is a true steady state that per- 
sists forever. The dashed blue line in Fig.[9]is a plot of the 
stress as a function of time for the full STZ model with 
a small initial perturbation to the effective temperature 
field, but the initial conditions are such that the defor- 
mation relaxes towards homogeneous flow. The simple 
ODE model stress solution, shown in magenta, lies on 
top of the PDE stress solution - the macroscopic stress 
response is the same for both models. The colored sym- 
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bols correspond to plots in Figs. fTOT a) and [TOfb) . 

Homogeneous deformation is characterized by the dis- 
sipation of perturbations to the effective temperature 
field. The effective temperature as a function of position 
is shown in Fig. llOf a). and each colored line represents 
the state of the system a different time. A small ini- 
tial perturbation to the effective temperature dissipates 
as a function of time, although the average value of the 
effective temperature increases as plastic work is dissi- 
pated. The effective temperature never varies more than 
5% from its average value. 

A similar plot for the plastic strain rate is shown in 
Fig. [TUf b). The plastic strain rate is zero during the 
elastic response, and although the perturbation to the 
effective temperature generates a small perturbation to 
the strain rate at the onset of plastic deformation, the 
strain rate relaxes towards a homogeneous state. The 
maximum plastic strain rate is remains within 20% of its 
average value. 
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FIG. 9: (color online) Deviatoric stress s vs. strain calculated 
by numerically integrating the STZ equations of motion with 
initial conditions Xim = 0.0674, q = 1.015 x 10"''. The dashed 
(blue) curve represents the solution to the perturbed system, 
while the dash-dotted (magenta) curve represents a homoge- 
neous solution where the effective temperature is constrained 
to be constant inside the material. In this plot the two curves 
are indistinguishable. The colored symbols correspond to the 
plots shown in Fig. 1101 Because Xim is large, the system be- 
gins with a large number of plasticity carriers and therefore 
the stress peak is negligible. This system does not localize. 
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FIG. 10: (color online) Relaxation to homogeneous flow: (a) 
effective temperature and (b) normalized plastic strain rate 
7(t/)ro/g as a function of position, for a material with initial 
conditions J dyx{y,t = 0) = Xini ~ 0.20, and imposed strain 
rate q = 1.015 x 10~^. Different (colored) lines represent 
different times. A small initial perturbation to the effective 
temperature dissipates, although the average effective tem- 
perature increases. (The effective temperature scale is much 
smaller than Fig. [S^a)). Initially the deformation is purely 
elastic and the plastic deformation is zero, and at the onset of 
plastic deformation the average plastic strain rate increases 
rapidly. Although the plastic strain rate is perturbed at this 
point, the perturbation decays. 



C. Diffusion limited shear bands and failure 

A second slowly evolving state, called "diffusion lim- 
ited localization," occurs when the shear heating and dif- 
fusion terms in Eq. ([5]) balance. In this case the effective 
temperature field is far from its steady state value x &t 
all points in space, so that the factor (1 — x/x) is close 
to unity and the shear heating term sx/{soCo), which is 
of order one, balances the diffusion term a^. The balance 



is not perfect, and x is not exactly zero, so the band 
continues to diffuse slowly outward. 

This type of deformation is important because it sets 
the minimum length scale for shear bands 0{a) ~ 0.015, 
in numerical STZ solutions. Although a subset of initial 
conditions generates shear bands that become thinner 
than this length scale, the local strain rate in these bands 
becomes so large that the steady state effective temper- 
ature X approaches infinity. In other words, sx/ (sqCo) 
is too large inside the band, the diffusive flux can not 
balance it, resulting in a runaway heating process. This 
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failure is not an artifact of our numerical methods; it sig- 
nifies a break down in the STZ model that occurs when 
the effective temperature increases without bound. We 
associate this runaway process with a third state, the 
onset of material failure, because the solid-like STZ de- 
scription fails as the material liquefies. 

Although we can not track the thickness of the band 
below the grid resolution during this runaway process, we 
do observe that just prior to failure the effective temper- 
ature in these simulations is elevated significantly above 
its average in a region of thickness a. In other words, 
the diffusion length scale appears to be an upper bound 
on the size of the region where structural changes occur 
during these shear failure events. 

Excluding material failure, the thinnest shear bands 
possible in this model are diffusion limited. These types 
of shear band persist for long times in STZ simulations for 
earthquake faults, where the parameters are highly rate 
weakening and chosen to reproduce the rate dependence 
observed in granular fault gouge experiments [2al ■ Diffu- 
sion also limits the initial thickness of shear bands in sev- 
eral of the numerical simulations performed in this paper, 
although these shear bands continue to diffuse outward 
at larger strains. 

Diffusion limited shear bands can be identified by their 
narrow thickness, which is of order a, although the ex- 
act value varies with the stress overshoot and specific 
heat Cq. In this 2D model, we use the term "thickness" 
to refer to the extent of the shear band in the direc- 
tion orthogonal to the slip plane, which is similar to the 
meaning of this term in three-dimensional systems 2]. 
Although diffusion of potential energy (and presum ably 
effective temperature) has been seen in simulations [2l|, 
the length scale a associated with this diffusion constant 
is relatively unconstrained by simulations or experiments. 
A reasonable postulate is that a is on the same order as 
the radius of an STZ, or equivalently, a few particle radii. 
This suggests that diffusion limited shear bands are very 
narrow. 

The stress vs. strain curve for a material that develops 
a diffusion limited shear band is given by the dashed blue 
line in Fig. [51 The stress weakens very rapidly as the 
diffusion limited shear band forms. For comparison, the 
dash-dotted magenta curve in Fig.[2]is the stress response 
of the ODE rate and state model with the same average 
initial conditions. This illustrates that thin shear bands 
drastically change the macroscopic system response, and 
that this dynamic weakening is not captured by a single 
degree of freedom rate and state model. 

As discussed in Section IIIIl Figures [3] and [5] show the 
time evolution of a shear band which is initially diffusion- 
limited. 



D. Disorder limited shear band 

The fourth "disorder limited" localized state is less in- 
tuitively obvious, but occurs frequently in our numerical 



STZ solutions. Neglecting the diffusion term, the right 
hand side of Eq. ^ is proportional to the product of 
two factors, exp(— 1/x) and (1 — x/x)- The former is 
very close to zero whenever x is significantly less than 
XOj and the latter is zero when x = X- The disorder lim- 
ited state occurs exactly when the small-x condition is 
met outside the shear band and X = X inside the band, 
so that X in Eq. ([5]) is always small. However, it is never 
zero, so that the disorder limited shear bands are also 
transient solutions that diffuse outward with time. This 
type of shear band was first described in 0, and cap- 
tures features of shear bands observed in simulations by 
Shi, et al. [2l[. Figure [TT] is a plot of the shear stress s 
vs. strain for a system that develops a disorder-limited 
shear band. This plot is calculated by numerically inte- 
grating the STZ equations of motion with initial condi- 
tions Xini = 0.1042 and ^ = 8.7 x 10"''. The blue curve 
represents the solution to the perturbed system, while 
the magenta curve represents a homogeneous solution 
where the effective temperature is constant as a func- 
tion of position inside the material. The colored symbols 
correspond to the plots shown in Fig. [TH Although the 
localized system weakens slightly faster than the homo- 
geneous system, the effect is small and on this scale the 
two curves are indistinguishable. 

We show the effective temperature and strain rate 
fields for a numerical solution that exhibits a disorder 
limited shear band in Figs. [T^a) and [T^ b). The per- 
turbation to the effective temperature field grows very 
slowly at first, then more rapidly as x — > x, and finally 
the peak begins to diffuse slowly outward. Similarly, the 
normalized plastic strain rate begins at zero (blue), then 
rises quickly (green) and relaxes slightly (red). 
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FIG. 11: (color online) Shear stress s vs. strain calculated 
by numerically integrating the STZ equations. The dashed 
(blue) curve represents the solution to the STZ PDE, while 
the dash-dotted (magenta) curve represents a solution where 
the effective temperature is constant inside the material. The 
colored symbols correspond to the plots shown in Fig. 1121 
Although the localized system weakens slightly faster than 
the homogeneous system, the effect is small and on this scale 
the two curves are indistinguishable. 
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FIG. 12; Disorder limited shear band(color online) (a) 
Normalized plastic strain rate 7(i/)To /q and (b) effective tem- 
perature as a function of position (y), for a material with 
initial conditions Xim = 0.1042, and imposed strain rate 
q = 8.7 X 10~®. Different (colored) lines represent different 
times; cooler colors (blue) correspond to earlier times, while 
warmer colors (red) correspond to later times. The plastic 
strain rate in the band increases significantly (about 800 %), 
although much less than in the diffusion limited shear band. 
The thickness of this band at its peak is about 0.2, much 
larger than the thickness of a diffusion limited shear band. 

The thickness of disorder limited bands is not set by a 
simple internal length scale such as a. Instead, the thick- 
ness is determined dynamically by the externally imposed 
strain rate and the initial conditions. 

Assume for the moment that a single shear band forms 
in the material. This is explicitly enforced for the numer- 
ical integration in this paper because the initial hyper- 
bolic secant perturbation at y = leads to a single shear 
band at that position. In addition, a single shear band is 
observed in simulations [2ll | and numerical integration of 
the STZ model with random perturbations to the initial 
effective temperature 0] at low strain rates. 

Under this assumption, almost all of the deformation 
is accommodated in a band of thickness w. 



Using Eqs. ^ and ([27)) we derive the following relation- 
ship between the stress s, the thickness of the shear band 
w and the externally imposed strain rate q: 



q2L 



w 



2/(s) exp 



1 



x{q2L/w) 



(28) 



This is not a prediction for the thickness of the shear 
band, because the final stress, s is not specified. Unfor- 
tunately, we can not derive an approximate value for s 
because it depends on the entire history of deformation in 
the material. In addition, the final value of s is generally 
close to the yield stress, and /(s) is very sensitive to s in 
this regime. However, in the next section we will check 
to see if the shear bands in a given numerical simulation 
satisfy the criterion given by Eq. 



E. 



Deformation map at the time of maximum 
deformation rate 



We now determine which of these states occur and per- 
sist as a function of the initial conditions in the numeri- 
cally integrated solutions. First, we use the following cat- 
egories to characterize the deformation at the time tmax ■ 
homogeneous deformation, diffusion limited shear band, 
disorder limited shear band, or material failure. This 
categorization is somewhat arbitrary because none of the 
states are stationary; all perturbed states will eventually 
decay towards homogeneous flow. However, by identi- 
fying these different regions in phase space we hope to 
identify length scales and features that are observable in 
experiments. 

To determine if a shear band thickness is consistent 
with disorder limited deformation, we rearrange Eq. (|28|). 

inserting Wqmax and Sqmax- 



log 



q2L 



1 



0. (29) 



A shear band in a numerical solution is said to be 
"disorder limited" if Eq. (|29)) is satisfied to within 8 %, 
(i.e., the left hand side equals ± 0.08.) Similarly, a 
shear band is "diffusion limited" if its thickness is ap- 
proximately equal to the diffusion length scale a, (i.e. 
< Wqmax < 0.03), and homogeneous if the maximum 
Gini coefficient is less than 0.5. These cutoffs are chosen 
to ensure that deformation regions are non-overlapping, 
which is a strong constraint. Finally, a material is said 
to fail if X approaches infinity during the course of inte- 
gration. Transition regions are expected when an inho- 
mogeneous flow does not flt into one of these categories. 
Figure[T3]is a deformation map that indicates where each 
of these criteria are satisfied. 



F. Deformation map at 20 % strain 



q = TqVo/L ~ To {w/L) jband- 



The same criteria for deformation categories that were 
(27) used at tqmax in Fig. [T3] can also be used to categorize 
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FIG. 13: (color online) Deformation map that uses the thick- 
ness shown in Fig. [7] to determine if the deformation at time 
tqmax is diffusion or disorder limited localization. Diffusion 
limited shear bands (very dark gray /dark red) and failure 
(black) occur where < it; < 0.035, and disorder limited shear 
bands (medium gray/orange) occur where the left-hand side 
of Eq. (I29|l is less than 0.7. The very light gray (light yellow) 
region indicates homogeneous flow. Because this is a snapshot 
of the the deformation types at t = tqmax, for some initial 
conditions the system is transitioning between two types of 
flows. The red region represents a transition regime between 
diffusion limited and disorder limited shear bands, while dark 
yellow represents a transition between disorder limited shear 
bands and homogeneous flow. Blue outline boxes indicate 
initial conditions detailed in Figs. 1101 and 1121 



FIG. 14: (color online) Deformation map at 20 % strain, 
that uses the thicknesss shown in Fig. [8] to determine the 
type of deformation. Failure (black) occurs where x ^ oo 
during a numerical simulation, and disorder limited shear 
bands (medium gray/orange) occur where the left-hand side 
of Eq. (|29|) is less than 0.7. The very light gray (light yel- 
low) region indicates homogeneous flow. At t = 0.2, the 
dark gray (red) transition region in Fig. [T3] has disappeared 
- the shear bands have widened to become disorder limited 
shear bands. In addition, the diffusion limited shear bands 
have also widened to become disorder limited, and some of 
the shear bands which were disorder limited in Fig. [13] have 
transitioned towards homogeneous flow (light gray /dark yel- 
low). Blue outline boxes indicate initial conditions detailed 
in Figs.OIISl andfT^ 



shear bands at 20% strain. Inhomogeneous flows which 
were in a "transition regime" at tqmax should broaden 
quickly towards one of the slowly varying categories. This 
is seen in Fig. [TH which shows deformation categories 
at 20 % strain. The transition region between disorder 
limited shear bands and failure has shrunk considerably. 
However, all of the diffusion limited flows have transi- 
tioned to disorder limited bands, and some disorder lim- 
ited shear bands have transitioned towards homogeneous 
flows. This highlights the fact that shear band thickness 
and deformation type depend significantly on the amount 
of strain. Most experiments are limited to small strains 
and therefore can not see this time evolution. 

This analysis of the STZ model shows that for a large 
range of initial conditions, shear bands are a robust fea- 
ture that persist for long times. STZ theory indicates 
that shear band thickness evolves slowly over very large 
strains and suggests that the thickness is determined dy- 
namically by the initial and boundary conditions. STZ 
theory predicts that the minimum thickness of the bands 
is set by an effective temperature diffusion parameter a, 
but a continuum of other thicknesses is also possible and 
dynamically determined. 

While we predict that the types of deformation 



mapped in Figs. [12] and [TJ] will occur in a wide range 
of amorphous solids, the exact location of boundaries be- 
tween types and the longevity of each type are likely ma- 
terial dependent. These deformation maps depend on the 
definition of the steady state effective temperature x{q). 

All deformation maps in this work apply specifically 
to the model glass simulated by Haxton and Liu because 
we used a function x{q) that fits their data for simulated 
repulsive disks; different materials may have slightly dif- 
ferent steady state effective temperatures, although it is 
also possible that x{q) is universal. The steady state ef- 
fective temperature can be measured in simulations by 
comparing the fluctuations and linear response of an ob- 
servable such as the pressure; this input is all that is 
needed to generate a deformation map for a new mate- 
rial using STZ theory. 

In the analysis, we have assumed a single shear band. 
Experiments and simulations of bulk metallic glasses 
show that the material develops multiple shear bands at 
higher strain rates [ll|, [2l[ . Developing a model for the 
number and spacing between shear bands is beyond the 
scope of this paper, but the STZ model should provide 
an excellent starting place for these analyses. 
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V. CONCLUSIONS 

We have analyzed the stabiHty of the STZ model with 
a strain rate dependent effective temperature, and found 
that the details of the rate dependence specify the steady 
state stability of homogeneous flows. Most simulated 
glasses exhibit rate strengthening, where the steady state 
stress increases as a function of strain rate, and we have 
shown that these materials are stable with respect to per- 
turbations in steady state. In contrast, rate weakening 
materials, such as granular fault gouge, are unstable with 
respect to shear bands in steady state. 

Perhaps surprisingly, shear bands develop even in rate 
strengthening materials. They result from an instability 
that develops during a transient stress response, when 
a material is driven from rest or driven at a new ve- 
locity. Although the perturbations are unstable only 
for small strains, the resulting inhomogeneous effective 
temperature profiles xiu) ^tre nearly stationary states of 
the model equations of motion and therefore these shear 
bands persist for long times. 

By including information about the rate dependence 
of the steady state effective temperature, xil): show 
that the STZ model generates a deformation map that 
includes homogeneous deformation, thick "disorder lim- 
ited" shear bands, thin "diffusion limited" shear bands, 
as well as the onset of material failure. 

The shear bands that emerge spontaneously in the STZ 
model capture several important features seen in simu- 
lations and experiments. First, the STZ model predicts 
that shear band formation coincides with stress relax- 
ation after the initial stress overshoot in start-up flows. 
In cases where the material does not fail, the model pre- 
dicts that shear bands gradually broaden over large ( 
> 20%) strains. 

For a fixed initial effective temperature, the STZ model 
predicts that the shear bands become thinner and that 
their internal structure becomes more disordered as the 
strain rate increases. This is similar to the "ductile to 
brittle" transition seen in amorphous materials as a func- 
tion of the strain rate [ll'| . At lower strain rates the ma- 
terial deforms nearly homogeneously and appears ductile, 
but at higher strain rates all the deformation is localized 
in a thin shear band or mode II crack. 

At very high strain rates and low initial effective tem- 
peratures, the effective temperature approaches infinity 
at the center of the band during the transient response 
and the system "melts" . The liquefied region is very thin 
and failure occurs near the maximum stress overshoot, 
which is consistent with material failure via shear band- 
ing seen at high strain rates in bulk metallic glasses. 

The model predicts that for materials that fail via 
this shear banding mechanism, the apparent shear band 
thickness should be at most the diffusion length scale a. 
Although a has not been measured experimentally, a rea- 
sonable assumption is that it is on the order of an STZ ra- 
dius. In bulk metallic glasses this scale should be at most 
30 atomic radii, on the order of 10 nm T|, which is much 



smaller than the thermal diffusion length scale (100-240 
nm 0), and this could explain the shear band thickness 
measured in these materials. While the STZ radius has 
not been estimated in granular fault gouge, this mech- 
anism could provide an explanation for the scale of the 
prominent fracture surface, which is orders of magnitude 
smaller than other length scales in earthquake faults. 

We have also shown that these localization dynamics 
can not be captured by a single degree of freedom rate 
and state friction law, and that analyzing steady state 
model dynamics can often be misleading. This is be- 
cause the structural degrees of freedom parametrized by 
x{y) evolve much more slowly than the stress dynamics, 
so that the microstructure continues to evolve although 
the stress appears to have reached a steady state. This 
insight is particularly important for materials that de- 
velop highly localized shear bands, as the friction law 
based on homogeneous dynamics is vastly different from 
one that accounts for transient shear band development. 
We suggest that localization may play a role in dynamic 
weakening seen at high shear speeds in granular mate- 
rials, and that the STZ PDE model generates a useful 
friction law in this case. 

While this is an exciting starting point for studying 
deformation and failure for amorphous materials, many 
fundamental questions remain. We discuss a few of them 
below. 

What is x{q) for various amorphous materials? 
Throughout this paper, we have used a fit to data gen- 
erated by Haxton and Liu as the definition for xio)- 
Haxton and Liu simulate a 2D amorphous packing of har- 
monically repulsive discs at thermal temperatures above 
and below the glass transition temperature, and used 
FDT to extract an effective temperature at each ther- 
mal temperature and strain rate. To our knowledge, this 
is the only such data set. It would be very interesting 
to use FDT to extract effective temperatures from sim- 
ulations of other types of amorphous packings, such as 
the Lennard Jones glass studied by Shi et al. [21[ , foams, 
amorphous silicon, or bulk metallic glasses. Is x{q) sim- 
ilar for all of these materials? Is the effective glass tran- 
sition temperature, xo, universal? 

One possibility is that the transition from glassy be- 
havior to simply activated behavior should occur when 
(7=1, (i.e., when the strain rate is the same as the inter- 
nal rate 1/to). However, it is also possible that in compli- 
cated materials like bulk metallic glasses, the transition 
occurs at slower rates than I/tq, since the STZs in these 
systems are large, multi-component regions that likely 
evolve more slowly than the phonon frequency. 

Are there other ways to measure x{q), such as 
looking at the behavior of a tracer harmonic oscillator 
inside a simulation box? Is it possible to define the 
effective temperature by quantifying the change in 
configurational entropy as a function of the potential 
energy? Numerical results from Ono, et al. suggest 
that this type of calculation is possible, but they were 
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not able to sample enough low probability states to 
state conclusively that the FDT and entropic definitions 
generate the same effective temperature. These are im- 
portant questions because it is very difficult to measure 
fluctuations in position or stress precisely enough in ex- 
periments to extract an effective temperature using FDT. 

What are the effects of geometry and boundary condi- 
tions ? We have so far restricted ourselves to the simplest 
possible shear geometry and symmetric, no conduction 
boundary conditions on the effective temperature. The 
boundary conditions on the effective temperature help 
determine the location of shear bands within the mate- 
rial as well as the steady states of the system. In many 
experiments and in some earthquake faults, shear bands 
tend to localize along the boundary ^26,] . Why does this 
occur? 

Different geometries can be modeled in STZ theory 
by adjusting the boundary conditions on the effective 
temperature. For example, a crystalline solid boundary 
might impose a constant, more ordered boundary condi- 
tion on the effective temperature, while a rough, jammed 
solid might do the opposite. It would be very interest- 
ing to investigate the effects of these conditions on shear 
band evolution. 

In addition, many engineering materials are tested 
under tension and compression, or a "notch" is placed 
on the surface of the material. In these cases there is a 
free boundary which can deform, leading to a coupling 
between deformation and stress. The necking instability 
has been investigated using earlier STZ models [23| 
- it would be interesting to repeat this analysis with 
our improved understanding of the coupling between 
structure and deformation. 

What is the connection between stick slip instabilities 
and shear banding? 

In Sectionini we showed that the steady-state stability 
of homogeneous flow was dependent on whether the 
material was rate strengthening or rate weakening - 
rate weakening materials were unstable with respect 
to shear bands. Interestingly, in rate-and-state (ODE) 
friction models, stick-slip instabilities can only occur 
when the system is rate weakening 28[. In addition, 
formation of a shear band coincides with what looks 
like a slip event in the macroscopic stress-strain curve. 
Can a shear band be understood as the PDE analogue 
of a slip event in a single degree of freedom ODE? Is is 
possible for the transient shear bands seen in ostensibly 
rate strengthening materials to generate stick-slip like 
behavior? Preliminary numerical solutions suggest that 
the transient shear bands can not generate stick slip 
behavior, but more work is needed on this avenue of 
research. 

Is the localized state weaker in absolute terms than the 
homogeneous state? Does this matter? Recently, re- 
searchers studying friction in fault gouge have found that 



the shear stress supported by the gouge weakens rapidly 
at high driving rates |[29| . Strain localization in the STZ 
model generically leads to a rapid decrease in the shear 
stress, and has been suggested as a mechanism for this 
experimental observation [25]. However, in most numeri- 
cal solutions to the STZ equations (e.g. Fig. [2), the final 
stress state in the localized system is equal to or higher 
than the stress in the homogeneous system, except in the 
special case where the initial effective temperature per- 
turbation is a step function (25| . In absolute terms, the 
localized system is stronger (or at least no weaker) than 
the homogeneous system, which is counter-intuitive. 

There are several ways to reconcile this information 
with intuition. First, we note that the rate at which 
the localized system weakens is much more rapid than 
the homogeneous system. For dynamic phenomena, such 
as stick-slip instabilities and stop-start experiments, the 
weakening rate and the total stress drop help determine 
the dynamic response. Is the rapid weakening seen in 
systems that localize large enough to cause stick-slip? 
Another possibility is that many of these systems attain 
strain rates at which the STZ solid-like description breaks 
down. Although we do not explicitly model this here, 
it seems likely that the liquid-like material in the band 
possess a vastly reduced strength compared to the solid 
outside the band. 
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APPENDIX A: STZ MODEL DETAILS 

A mean field theory for shear transformation zones has 
been developed in a series of papers and 
we use this theory as a general model for a wide range of 
amorphous solids. 

In analyzing the dynamics of shear transformation 
zones, we develop equations of motion for five internal 
variables: the deviatoric stress s, the pressure p, the den- 
sity of STZs oriented parallel and perpendicular to the 
principal stress directions n± , and the effective tempera- 
ture, X- In a simple shear geometry at low temperatures, 
the model can be further simplified so that the state of 
the system is entirely specified by s and x alone. The 
following sections review the STZ equations and specify 
the parameters and simplifications used in this paper. 
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1. Overview of equations of motion 

In the slowly sheared materials we are modeling, the 
speed of sound in the material is very fast compared to 
the rate of plastic deformation. In this case the stress 
gradients equilibrate very quickly, and we take the zero 
density limit of the momentum conservation equations. 
This results in static elastic equations for the stress: 



= 0. 



(Al) 



The rate of deformation tensor is the sum of elastic and 
plastic parts: 

1 / dvi , dvi 
2 



D. 



total 



dxj ~^ dxi 



where V/'Dt is the material or co-rotational derivative. 
To simplify notation, the deviatoric stress has been 
nondimensionalized by an effective shear modulus Sy that 
specifies the stiffness of the STZs. The stress scale Sy also 
characterizes the stress at which the material begins to 
deform plastically. This yield stress is distinct from the 
maximum stress attained, Sm, and the steady state flow 
stress, s/, both of which are sometimes also referred to 
as the yield stress in the literature. 

The plastic rate of deformation tensor can be written 
in terms of dynamical variables from STZ theory. We 
postulate that under shear stress, each STZ deforms to 
accommodate a certain amount of shear strain, and can- 
not deform further in the same direction. This is modeled 
by requiring that each STZ be in one of two states: ori- 
ented along the principal stress axis in the direction of 
applied shear, which we will denote , or in the per- 
pendicular direction, "— " . 

Under applied strain, the STZ will flip in the direction 
of strain, from "— " to . Under shear stress in the 
opposite direction, the STZs can revert to their original 
configurations, which corresponds to a flip from to 
"— " . We assume that the STZ density is small and each 
STZ interacts with other STZs through continuum fields 
such as the stress. Therefore the rearrangements or flips 
occur at a rate R{s)/to, which depends on the stress and 
a characteristic attempt frequency I/tq. 

Because each STZ can flip at most once in the direction 
of applied strain, STZs must be created and annihilated 
to sustain plastic flow. Based on these considerations, 
the number density of STZs in each direction, n±, obeys 
the following differential equation 



n± 



Ton± = R{±s)n^ - R{Ts)n± + T (^^e'^^^ 

(A3) 

where R{±s)/to is the rate of switching per STZ as a 
function of stress, T is the rate at which energy is dissi- 
pated per STZ, and e~^^^ is the steady state density 
of STZs in equilibrium. 



The plastic rate of deformation tensor is given by the 
rate at which STZs flip: 



nooTo 



(i?(s)n_ - R{~s)n+) , 



(A4) 



where eo is a strain increment of order unity and rioo is 
a density roughly equal to the inverse of the volume per 
particle. 

The first two terms in Eq. (jASP correspond to STZs 
switching from to "— " states and vice-versa, while 
the last term enforces detailed balance: STZs are created 
at a rate proportional to noo^~^^^ and annihilated at a 
rate proportional to their density. The creation rate is 
proportional to the probability of a configurational fluc- 
tuation that corresponds to an STZ. As discussed in the 
introduction, this probability is exp[— 1/x], where x is 
an internal state variable that characterizes the conflgu- 
rational disorder. To close the system of equations, the 
model requires an equation of motion for y. 

Ono et al. jT3| and Haxton and Liu Q show that a 
driven amorphous system possesses a well-defined steady 
state effective temperature, x, at each value of the im- 
posed strain rate. In these simulations, a thermostat 
ensures homogeneous deformation within the glass and 
the particles are sheared for long periods of time before 
the steady state measurement is taken. 

To model deformation in time varying systems, such as 
start-up flows, we have to estimate how the effective tem- 
perature changes in time. As detailed in p7| . we assume 
that the heat content in the configurational degrees of 
freedom is driven by two independent sources, mechan- 
ical work and thermal fiuctuations. Therefore we make 
the simplest assumption: the mechanical heat drives the 
effective temperature towards x according to the con- 
ventional linear law of heating. The rate of heat per unit 
volume that enters the configurational degrees of freedom 
is Q^=Teff {dSc/dt)^^^f^. 

In addition, we postulate that the heat produced by 
thermal fluctuations, Qt = T^f f {dSc/dt)^f^^^^, drives 
the effective temperature towards the thermal equilib- 
rium bath temperature 9 according to the linear law of 
cooling. The resulting equation of motion for x is: 



1 



+ Teff 



Teff 

dSc 
~dt 



dSc 
~dt 



the 



1 



T 



where Ceff is a specific heat, Tz = Ez/ks is the STZ 
formation energy in temperature units, and the last term 
represents diffusion of effective temperature. 

Because the effective temperature governs the config- 
urational degrees of freedom, only configurational rear- 
rangements, i.e. plastic events, permit diffusion of the 
effective temperature. This suggests that the diffusivity 
should vary with strain rate, so that D = a^|7p;|, where 
a is a length scale that corresponds to the radius of an 
STZ. 
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2. Simplifying assumptions 



ergy dissipation rate, Q, is proportional to F: 



Pechenik [331 generalized Eqs. (|X3)l and (jM)) to the 
case where the principal axes of the STZ orientation ten- 
sor riij are not aligned with principal axes of the stress 
tensor Sij . These generalized equations can be written in 
terms of two new variables, A and m, which appear often 
in literature on STZs: 



l{j — Tl{j /tIqq , 



(A6) 
(A7) 



where ritot is the tensorial generalization of (rt+ + 
and Tiij is the tensorial generalization of («+ — rt_). The 
scalar A is the total density of zones in a sample, while 
the tensor m corresponds to the STZ orientational bias. 

In this paper we focus on materials in a 2D simple 
shear geometry, so that the diagonal terms in the devi- 
atoric stress tensor {s^x, Syy) and STZ orientational bias 
{rrixxj'myy) are significantly smaller than off-diagonal 
terms and can be neglected. Let s = s 



Syx and 



m = TTixy = ruyx- In this geometry the pressure p does 



not change with time and Vp/Vt in Eq. (jA2p is zero. 

As noted in ^S^j the density of STZs, eoA, is necessarily 
small. In a simple shear geometry, the equations of mo- 
tion for the stress s and the effective temperature x each 
contain this factor in their numerators, and they equili- 
brate very slowly compared to m and A. Therefore we 
replace A and m by their steady state values. Combining 
Eqs. ((X3|) . (|X6|l . and (|T7)) . we find that the steady state 
value of A is exp[— 1/x], and that m exchanges between 
two steady states (elastic vs. plastic deformation) near 
s = Sy. Below the yield stress the deformation is almost 
entirely elastic because all the existing STZs are already 
flipped in the direction of stress. Above the yield stress 
STZs are continuously created and annihilated to sustain 
plastic flow. Details can be found elsewhere (l7| . 

We also make the simplifying assumption that the ma- 
terial is below the thermal glass transition temperature 
To, so that particle rearrangements are not activated by 
thermal fluctuations. In this case the thermal entropy 
contribution to the effective temperature equation of mo- 
tion, (rf'5'c/rfi)t;jerm ^Iq. (|A5|) . is zero. This is always 
true, for example, in granular materials. In addition, we 
make the approximation that at very low temperatures 
the STZs do not flip in a direction opposite the direction 
of apphed stress: |s|) — 0. 

The functional form of F, the energy dissipated per 
STZ that appears in Eq. (|A3p . is considerably simplified 
for T < Tq. Under the assumption that an STZ does 
not flip in a direction opposite to the direction of applied 
stress, no energy can be stored in the plastic degrees of 
freedom. This means that the plastic work is equal to 
the energy dissipated: 



D^'s- - 



7S = Q- 



(A8) 



Following Pechenik [33|, we postulate that the total en- 



Q = so-AF, 

To 



(A9) 



where so is a stress scale we return to below. 

Combining Eas. lA4[|A8[ and lAQI results in the following 
expression for F: 



F(s) 



(AlO) 



where the function /(s), which also appears in Eq. ([2]), 
is defined as follows: 



/(^) 



eo- 



7^(s) 



1 — m(s) 



(All) 



In simple shear below the thermal glass transition, the 
two steady states of the equation of motion for the STZ 
bias are also very simple: 



, , If for s < Sy 

Ms) < , r ^ 

I So/ s for s > Sy, 



(A12) 



and So = Sy is the yield stress - the stress at which the 
two stability branches intersect. 

We now turn to the parameters in Eq. (|A5p . In [l3| the 
rate at which configurational entropy is being produced 
by mechanical deformation, [dSd dt)mech, is assumed to 
be proportional to the product of the STZ density and 
the STZ creation rate: 



dSc 
~dt 



Vl To' 



AF, 



(A13) 



where $7 is a volume per molecule and ksvz is the entropy 
per STZ [I^. 

Putting everything together, we arrive at an equation 
of motion for the effective temperature: 



dx ^ 2 sx 
CQSoq- 



xiq) 



where 7 is strain, q is the imposed strain rate times the 
STZ time scale, (Vb/i)To and cq = Ceff^l [kBVz)- In- 
serting Eq. into the equation for the rate of deforma- 
tion tensor, Eq. (|A2p . and integrating across the width 
of the material (in the y-direction) results in a second 
equation for the stress dynamics: 



(A15) 



where /i» is the ratio of the elastic modulus to the yield 
stress, and A is the spatial average of the STZ density 
A = exp(-l/x). 

Determining an exact function R{s) from first prin- 
ciples is a difficult many body problem. However, we 
do know how R{s) behaves in the limits of very small 
and very large stresses, and we choose a function that 
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smoothly interpolates between these two regimes. R{s) 
exhibits Eyring-like behavior far below the yield stress, 
and power law behavior above the yield stress: 



The first factor on the right-hand side of Eq. (jA16p is the 
Eyring rate in a form similar to that used in [16i] . where 



the exponential function of s/p, causes the rate to sat- 
urate at large s. Here, is the height of the Eyring 
activation barrier in units of temperature. The expo- 
nent n in Eq. (jA16[) specifies the large stress power law 
behavior; possible values are discussed in Q. Analysis 
of simulation and experimental data suggest that ri = 1 
is valid for bulk metallic glasses [l3l but in this paper 
we use n — 1/2, which is relevant for purely repulsive 
harmonic disks (8| . 
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